Large eddy simulation of the bubbly ocean: New insights on subsurface bubble distribution and bubble-mediated gas transfer

[1] The evolution of bubbles in a turbulent oceanic boundary layer is simulated using a multi-size multi-gas bubble model coupled with a Large Eddy Simulation model. Bubbles injected by breaking waves are brought into the boundary layer by episodic bubble plumes, and they form near-surface streaks in the convergence zone of Langmuir circulations. The equilibrium bubble distribution decays exponentially with depth and is a manifestation of intermittent bubble plumes whose bubble number density is at least one order of magnitude higher than the mean bubble number density. Bubble distribution in the injection zone is influenced by injection, turbulent transport, and dissolution. Bubble distribution below the injection zone is determined by the strength of turbulence and dissolution. For a given sea state, bubble e-folding depth increases linearly with friction velocity. Wave age is an additional governing parameter for bubble e-folding depth. The buoyancy of bubbles weakens both Langmuir circulations and near-surface turbulent kinetic energy dissipation. The buoyancy effect increases with wind speed. Gas flux through bubbles depends on both wind speed and wave age. For a given sea state, the bubble flux increases with wind speed to the fifth power.